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Abstract 

The vector meson mass is extracted from a large sample of partially quenched, 
two-flavor lattice QCD simulations. For the first time, discretisation, finite-volume 
and partial quenching artefacts are treated in a unified framework which is consis- 
tent with the low-energy behaviour of QCD. This analysis incorporates the leading 
infrared behaviour dictated by chiral effective field theory. As the two-pion decay 
channel cannot be described by a low-energy expansion alone, a highly-constrained 
model for the decay channel of the rho-meson is introduced. The latter is essen- 
tial for extrapolating lattice results from the quark-mass regime where the rho is 
observed to be a physical bound state. 



Recent developments in lattice QCD have enabled the first large-scale simu- 
lation of chiral, dynamical fermions [1]. While this accomplishment is a sig- 
nificant milestone in the progress towards an accurate description of physical 
QCD, the high demand on computing resources restricts practical calcula- 
tions to an unphysical domain of simulation parameters. In particular, lat- 
tice QCD involves a discretised space-time of finite spatial extent, with input 
quark masses much larger than those in Nature. Each of these approxima- 
tions requires special attention in the extraction of physical observables from 
Monte Carlo simulations. In this Letter we analyse a very large set of partially 
quenched data for the mass of the p meson. We show that a systematic anal- 
ysis of this data enables us to remove the effects of partial quenching and to 
take both the continuum and infinite volume limits. The resulting data lies on 
a single, well defined curve which extrapolates to a value within ~1% of the 
physical p mass. The contrast between the raw lattice data (note the scatter 
in Fig. 2) and the corrected data shown in Fig. 3 is striking. 
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Spontaneous chiral symmetry breaking in QCD dictates that, in the vicin- 
ity of the chiral limit, hadronic observables exhibit nonanalytic dependence 
on the quark mass [2]. This feature places tight constraints on the form of 
chiral extrapolations if they are to be consistent with the properties of low- 
energy QCD [3,4]. The most natural solution to this problem is to use effective 
field theory (EFT) to describe the quark-mass dependence of hadron proper- 
ties. Considering a benchmark quantity, such as the nucleon mass, there is 
substantial phenomenological information on the quark-mass expansion near 
the chiral limit [5] . In the context of lattice simulations, where quark masses 
are significantly far from the chiral limit, the expansion is acutely sensitive 
to higher-order corrections. Fortunately, such issues can be alleviated by re- 
formulating the EFT in the framework of finite-range regularisation (FRR) 
[6] ■ - with demonstrated success in the efficient extrapolation [7] of lattice 
calculations of the nucleon mass [8] . 

Provided simulations are performed on a suitably large box, finite-volume 
corrections will be exponentially suppressed. Nevertheless, these leading cor- 
rections can be described by the same low-energy effective field theory used to 
understand the quark- mass variation [9]. Finite- volume corrections are domi- 
nated by the suppression of the infrared component of chiral loop diagrams - 
as observed in a recent study of volume dependence in lattice QCD [10] (using 
the quark-mass dependence described in Ref. [11]). Corrections of this type 
have previously been incorporated in quark-mass extrapolations [12,13,14]. 
They are essential in the case of the p-wave decay channels, such as p — > im 
[12] and A — > Nir [13]. The modifications to EFT on a finite volume have 
been investigated for a range of observables — e.g., see Refs. [15,16,17,18,19]. 

Removal of discretisation artefacts from simulation results also represents an 
important step in the systematic extraction of continuum QCD physics. From 
a technical point of view, a great deal of effort has gone into action improve- 
ment [20] such that near-continuum results can be obtained at finite lattice 
spacing [21,22,23]. Residual discrepancies from the continuum can be incor- 
porated as perturbative corrections in EFT [24,25,26,27], thereby providing a 
systematic approach to the continuum. 

The generation of gauge field configurations with dynamical sea quarks is the 
most computationally demanding component of the calculation of standard 
observables. By comparison, the matrix inversion required to obtain quark 
propagators is relatively efficient. This enables the calculation of quark prop- 
agators over a range of quark masses for a fixed gauge field ensemble. Such 
calculations are referred to as partially quenched QCD (pQQCD), where the 
valence quark masses no longer match those simulated in the sea. Although 
an unphysical approximation, the connection to the physical theory in EFT 
has been demonstrated [28]. Most importantly, the partially quenched EFT 
does not require any new, unphysical parameters. 
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While EFT provides a systematic framework for the analysis of lattice results, 
the present analysis of the p meson requires one to go beyond the scope of 
EFT. Near the chiral limit, the p decays to two energetic pions. The pions con- 
tributing to the imaginary part of the p mass cannot be considered soft, and 
therefore cannot be systematically incorporated into a low-energy counting 
scheme [29,30]. Because almost all the lattice simulation points in this analy- 
sis lie in the region m n > m p /2, it is evident that in the extrapolation to the 
chiral regime will encounter a threshold effect where the decay channel opens. 
To incorporate this physical threshold, we model the p — > -kit self-energy dia- 
gram constrained to reproduce the observed width at the physical pion mass. 
Including this contribution also provides a model of the finite volume correc- 
tions arising from the infrared component of the loop integral. In particular, 
we can also describe the lattice results in the region m, n < m p /2, where the 
decay channel still has higher energy because of momentum discretisation. 

In this Letter we present a global analysis of a very large set of sophisticated 
lattice simulation results. The aim is to produce an accurate determination of 
the (real part of the) p-meson mass in 2-flavour QCD, with the required input 
of the experimental width. For a complete account of the analysis procedure, 
the reader is referred to Ref. [31]. 

Partially-quenched lattice simulations are characterised by the distinction be- 
tween the masses of quarks coupling to external sources and those associated 
with vacuum polarisation of the gauge field. Subsequently, the construction of 
effective field theories based on such simulations necessarily distinguishes the 
valence- and sea- quark composition of hadronic states. While the external 
legs of any n-point hadronic correlator are constructed of valence quarks only, 
internal (hadron) loop diagrams may comprise any mixture of sea and valence 
quarks. The mass of a particular vector meson state is therefore described by 

M? jk = M(a,m i ;m j ,mk), (1) 

where the first two parameters of M(a,mi]mj,mk) specify the gauge field 
ensemble with lattice spacing a and sea-quark mass m sca = m^. The masses 
of the valence quarks are given by rrij and m^. Similarly, mfj k is the mass of 
a pseudoscalar meson of equivalent quark composition. 

The global parameterisation of the vector meson mass, dependent on quark 
masses, lattice spacing and physical volume, is written as 

M£/ = (ao + X 2 a 2 + a 2 m%? + a 4 m^) 2 + SjgJ(L) . (2) 

The total loop correction to the M?- ■ meson, on a finite box of physical length 
L = Na, is described by 

S™J(L) = L) + L) + Ej p (m? y , m^, m^, L) . (3) 
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Fig. 1. Left: Diagram providing the leading nonanalytic contribution to the chi- 
ral expansion of the p- meson mass (a) and its associated quark- flow (b). Middle: 
Two-pion contributions, (a), (c), to the p- meson self-energy and their associated 
quark-flow, (b), (d), respectively. Right: The rj contribution (a) and its associated 
quark flow diagrams in pQQCD. While diagram (c) appears in quenched QCD, in 
pQQCD (or full QCD) it is complemented by an infinite series of terms, the first 
two of which are depicted in diagrams (d) and (e). 

The relevant loop corrections are depicted in Fig. 1. The corresponding loop 
integrals, in the L — > oo limit, are given by 
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Symbols are summarised as: = 93 MeV; physical p and 7r masses, \x p = 
0.770 GeV and ^ = 0.140 GeV; pmr coupling f p7TW = 6.028 [12]; = [i p g\\ 
g>2 (as introduced in Ref. [32]) is related to the uipn coupling, g Wjm = 16 GeV -1 
[33], by g 2 = gup-nfir/^ = 0.74; k = \k\; finite-range regularisation is imple- 
mented with a dipole form, u(k) = (1 + k 2 /A 2 )~ 2 , and the pirn coupling is 
preserved at the physical threshold, u nn (k) = u(k)u^ 1 (^J /i 2 /4 — p, 2 ). By de- 
manding this physical threshold the infrared behaviour is no longer constrained 
and hence is not controlled by EFT. Therefore, such a model is essential to 
extrapolate lattice results from the regime > m p /2. 
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Table 1 

The resulting global fit to the entire sample of lattice results (with A = 655 MeV). 



a 


x 2 


a 2 


«4 


X 2 /d.o.f. 


(GeV) 


(GeV-fm- 2 ) 


(GeV- 1 ) 


(GeV" 3 ) 




0.8321* 


-1.4012 


0.49411? 


-0.06111 


39 / 76 



The leading finite-volume corrections are trivially incorporated by replacing 
the continuum loop integrals in Eqs. (4,5,6) by a sum over discrete momenta 
[12,13,14] 

k 

where k = (2n/L)i for % G Z 3 . This modification produces an infrared sup- 
pression of the loop integrals, and is independent of the choice of ultra-violet 
regularisation [16]. 

The bracketed term in Eq. (2) describes the residual variation of the vector 
meson mass which is not contained in the one-loop Goldstone boson diagrams. 
The analytic variation of the quark-mass dependence is characterised by the 
continuum parameters, ctj. At finite lattice spacing, all terms at order a and 
a 2 must be treated consistently with the symmetry breaking patterns of the 
prescribed fermion action [24]. This can potentially lead to more singular 
chiral behaviour in the effective field theory at finite lattice spacing [26]. In 
this study, the leading lattice spacing corrections to the terms analytic in the 
quark mass are investigated. 

The lattice simulation data considered in this analysis come from a large 
sample of partially quenched simulation results from the CP-PACS Collabo- 
ration [8]. These simulations were performed using mean-field improved clover 
fermions at four different couplings, (5. For each value of the coupling, four 
different sea quark masses have been calculated, yielding a total of 16 inde- 
pendent gauge field ensembles. On each ensemble, the quark propagator has 
been evaluated for five values of the valence quark mass. For the vector mesons 
constructed of degenerate valence quarks there are a total of 80 "data" points. 
The lattice scale is set via the QCD Sommer scale r = 0.49 fm [34], enabling 
all of these points to be shown in physical units — as in Fig. 2. 

The form provided by Eq. (2) allows a universal fit to all these points with 
just four free parameters (plus regulator scale). This is therefore a highly con- 
strained fit to the large sample of simulation results. The best fit parameters 
are displayed in Table I. The x 2 indicates that Eq. (2) accurately describes 
this large quantity of data. The regulator mass, A, has also been optimised to 
produce a best fit to the data, namely A = 655 ± 35 MeV, with the bound de- 
termined statistically from the 1-a variation about the central fit. Studying the 
variation of the fit over this domain introduces a small additional uncertainty 
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Fig. 2. Partially-quenched vector meson masses plotted versus "degenerate" pion 
mass squared. Here = Mf- and m^ eg ^ = m^j. The simulation results are 

from CP-PACS [8]. 



to the extrapolated result which is listed below in the error estimate. 

Variations of Eq. (2) have also been investigated. Scaling corrections to the 
parameters 0:2 and 0:4 yield coefficients which are consistent with zero. Linear 
corrections in the lattice spacing, taking the form of a term X±a, are observed 
to be small and hence do not improve on the fit with X\ = 0. Extending to 
higher analytic order in the quark mass expansion, by a term a^m®, reduces 
the stability of the fit, indicating that the data are consistent with a 6 = 
- see Ref. [31] for a complete account of these effects. The systematic error 
quoted below covers the range found with all of these variations. 

The fit parameters shown in Table 1 allow one to shift the simulation results 
to the infinite- volume, continuum limit and to remove the effects of partial 
quenching — hence restoring unitarity in the quark masses. Complete details 
of the procedure are outlined in Ref. [31]. The results are displayed in Fig. 3, 
where we observe a remarkable result. The tremendous spread of data seen in 
Fig. 2 is dramatically reduced, with all 80 points now lying very accurately 
on a universal curve. 



The curve through Fig. 3 displays the determined variation of the p-meson 
mass with pion mass. This curve also presents an extrapolation to the physical 
point, allowing extraction of the physical p-meson mass 

M p = 778(4)^ 6 (8) MeV, (8) 

where the first error is statistical, the second is systematic and the third from 
the determination of A [31]. This result is in excellent agreement with the 



6 




Fig. 3. The same 80 lattice data points as in Fig. 2, after correction to restore 
the infinite- volume, continuum and quark-mass unitarity limits. The central curve 
displays the best-fit from the global analysis. The dashed curves show the bounds 
on the FRR scale, 0.620 < A < 0.690 GeV. 

experimentally observed mass. 

The systematic uncertainty arises from the choice of fitting function, as out- 
lined above, and also from the choice of finite-range regulator. In addition to 
the presented dipole form, the analysis has been repeated with monopole and 
Gaussian regulators. The central values of the extrapolated result with the 
monopole and Gaussian forms differs by +3 and —6 MeV, respectively, from 
the dipole result. Each regularisation scheme produces a different model of the 
p — > mi vertex. This suggests that the model-dependence of this contribution 
is small, once constrained to produce the correct width at the physical point. 

This extrapolation of the lattice results also offers an estimate of the vector 
meson mass in the chiral limit, M°. The central value of the preferred fit gives a 
value 775 MeV, with errors similar to those quoted above for the vector mass 
at the physical point. These errors are from extrapolation of lattice results 
only, whereas for phenomenological purposes the physical point provides a 
further constraint. We therefore report the correlated difference between the 
physical and chiral limit value 

M p - M° = 3.7(2)3(8) MeV . (9) 

Interestingly, the chiral limit value of M p is very similar to that of the physical 
value. This feature is observed in the reduction in slope of the extrapolation 
curve in Figure 3 as the chiral limit is approached. The underlying physics giv- 
ing rise to this reduced slope is the presence of substantial spectral strength 
in the low-energy two-pion channel below the rho-meson mass [35]. This sug- 
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gests a small sigma term for the p in comparison with the nucleon, where the 
curvature is enhanced near the chiral limit [7]. 

This analysis demonstrates the ability to treat all lattice artifacts within a 
unified framework. Both scaling violations and finite- volume discrepancies can 
be removed through the procedure outlined. The number of simulation points 
can be increased dramatically by including partially quenched results. This in 
turn permits a highly constrained fit to produce an accurate extrapolation to 
the physical point. With minimal input, namely the pmr and upn coupling 
constants, the real part of the p-meson mass has been accurately determined 
in two-flavour QCD 1 . The final result for the pion mass variation, as described 
by the universal curve in Fig. 3, sets a benchmark for the continuum, infinite- 
volume limit of the p meson in two-flavour QCD. 
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